Report for Model 9

— for both jointly and separately

2024-08-20

1. Joint Modeling

1.1 Trace plots for convergence check

The current MCMC setting is:

  • 300,000 iteration;
  • 200,000 burn-in;
  • 100 thinning.

1.2 Gelman and Rubin’s convergence check

## Potential scale reduction factors:
## 
##         Point est. Upper C.I.
## HDevsum      0.999       1.00
## LDevsum      0.999       1.00
## dh0          1.015       1.05
## dh1          1.015       1.05
## dl0          1.868       3.52
## dl1          1.718       3.22
## dl2          1.336       2.05
## dl3          1.096       1.23
## 
## Multivariate psrf
## 
## 1.78

1.3 ACF Plots

Here we plotted ACF plots for the following variables:

  • Total deviance;
  • Variables that didn’t pass the convergence check.

1.4 WAIC results

LevelH LevelL
DIC 1224.97003 22727.5825
DIC3 1161.88355 20648.2725
PWAIC 44.02095 996.1293
WAIC 1189.68711 21027.5778

2. Separate Modeling

2.1 Trace plots for convergence check

The current MCMC setting is:

  • 300,000 iteration;
  • 200,000 burn-in;
  • 100 thinning.

2.2 Gelman and Rubin’s convergence check

## Potential scale reduction factors:
## 
##         Point est. Upper C.I.
## HDevsum       1.01       1.03
## LDevsum       1.01       1.03
## dh0           1.60       2.91
## dh1           1.52       2.73
## dl0           1.37       2.49
## dl1           2.57       5.95
## dl2           1.02       1.05
## dl3           2.43       6.74
## 
## Multivariate psrf
## 
## 1.87

2.3 ACF Plots

Here we plotted ACF plots for the following variables:

  • Total deviance;
  • Variables that didn’t pass the convergence check.

2.4 WAIC results

LevelH LevelL
DIC 1448.01985 23921.990
DIC3 1292.89309 21168.527
PWAIC 94.59621 1536.806
WAIC 1360.96116 21915.208